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ABSTRACT 

Non-ionizing stellar continua are a potential source of photons for continuum 
pumping in the hydrogen Lyman transitions. In the environments where these 
transitions are optically thick, deexcitation occurs through higher series lines. 
As a result, the emitted flux in the affected lines has a fluorescent contribution 
in addition to the usual recombination one; in particular, Balmer emissivities 
are systematically enhanced above case B predictions. The effectiveness of such 
mechanism in H II regions and the adequacy of photoionization models as a tool 
to study it are the two main focuses of this work. 

We find that photoionization models of H II regions illuminated by low- 
resolution (X/8X < 1000) population synthesis models significantly overpredict 
the fluorescent contribution to the Balmer lines; the bias has typical values of 
the order of a few hundredths of a dex, with the exact figure depending on the 
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parameters of the specific model and the simulated aperture. Conversely, pho- 
toionization models in which the non-ionizing part of the continuum is omitted 
or is not transferred significantly underpredict the fluorescent contribution to the 
Balmer lines, producing a bias of similar amplitude in the opposite direction. 

Realistic estimations of the actual fluorescent fraction of the Balmer inten- 
sity require photoionization models in which the relevant portion of the stellar 
continuum is adequately represented, i.e. its resolution is high in the region 
of the Lyman lines. In this paper, we carry out such estimation and discuss 
the variations to be expected as the simulated observational setup and the stellar 
population's parameters are varied. In all the cases explored, we find that fluores- 
cent excitation provides a significant contribution to the total Balmer emissivity. 
We also show that differential fluorescent enhancement may produce line-of-sight 
differences in the Balmer decrement, mimicking interstellar extinction. Fluores- 
cent excitation emerges from our study as a small but important mechanism for 
the enhancement of Balmer lines. As such, we recommend to take it into account 
in the abundance analysis of photoionized regions, particularly in the case of 
high-precision applications such as the determination of primordial helium. 

Subject headings: atomic processes — ISM: abundances — H II regions — line: 
formation — stars: atmospheres 



1. Introduction 



The systematic study of gaseous nebulae starte d in the decade of 1930 's with a long 
series of papers by M enzel and collaborators (e.g., iMenzell 119371 ; iMenzel fc Bakerl 119371 ; 
Baker &: Menzell Il938l . among others). In these papers, the basis was laid for the quan- 
titative analysis of the hydrogen spectrum and the determination of physical conditions in 
gaseous nebulae. The hydrogen spectrum was recognized to form in the radiative cascade 
following recombination of H + atoms, and quantitative predictions for the line intensities 
based on this scenario were made. Two limiting cases were considered, which depend on the 
optical thickness of the nebula in the Lyman lines. In case A, the nebula is assumed optically 
thin in all recombination lines, including those of the Lyman series; in case B, the nebula 
is assumed thick to Lyman radiation, which is consequently scattered until conversion into 
higher series photons plus Lya or 2-7 continuum emission takes place. Case B was soon 
shown to be the most ap propriate description for most nebulae , both on observational and 
theoretical grounds (e.g., 



Baker fc Menzell Il938l : |Pottaschlll960l ). 



A third scenario, named case C, was also initially introduced to represent the case of 
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optically thin nebulae with stellar continuum pumping as an additional mechanism for the 
excitation of Lyct. Contrarily to cases A and B, howev er, case C has barely had any follo w-up 



in the literatur e, with very sp oradic exceptions (e.g.. IChamberlainlll953l ; lFerlandlll999l ). As 



pointed out by iFerland a reason for this scarce popularity is the a priori assumption 



of the relative lack of photons at the relevant frequencies due to the absorption features 
of stellar atmospheres, which would make the mechanism inefficient and imply that case C 
might be a merely academical scenario with no actual physical applications. (The initial 
calculation by Baker & Menzel did not suffer from such limitation, since the stellar source 
was approximated by a blackbody.) 

A further reason for the oblivion of case C is that the scenario it modifies (case A) is 
inadequate to represent most detectable objects. In other words, the lack of interest toward 
case C might not have to do with the correction it applies to case A, but rather with case A 
itself; after all, that nebulae are irradiated by stellar light is beyond doubt. If this is the case, 
one might wonder why not considering the effect of stellar irradiation as a correction to case 
B, rather than case A. Since case C is the irradiated counterpart to case A, it seems natural 
to christen "case D" the irradiated counterpart to case B (Tabled]). Note that cases A and 
B depend solely on the physical conditions in the gas, while cases C and D also depend on 
the illuminating source, so they are not cases in the same 'local' sense of cases A and B. 

Case D has never been taken, to the best of our knowledge, into consideration. Although 
this might seem somewhat surprising at first glance, we believe that there are reasons for it; 
the arguments developed in this paper will help to understand why, so we defer a discussion 
of this point to the last section. 

In Section [21 we will argue that case D is a better representation of real nebulae than case 
B, the difference between the two scenarios being mostly small but substantial by modern 
standards of accuracy. Curiously enough, case D is implicitly assumed in many (but not 
all) photoionization models, although in a biased way; the bias can only be eliminated by 
feeding into the models high-resolution near-UV spectra (Section [3]). We will also show 
that the differential fluorescent contribution to the various Balmer lines can also affect the 



Table 1: Nomenclature for different asymptotic behaviors ("cases") of the hydrogen spectrum 
in gaseous nebulae. Column describers specify the optical depth in the Lyman lines, row 
describers specify whether irradiation by stellar continuum is included or not. 





thin 


thick 


Non-irradiated 


A 


B 


Irradiated 


C 


D 
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intrinsic Balmer decrement, pushing it toward either the blue or the red (Section Hj). The 
resulting spatial distribution of the Balmer decrement depends on the spatial distribution 
of stars and the gas around them. A brief discussion of the dependence of fluorescence 
enhancement on the nebular and the stellar population's parameters is given in Section 
and on details of the assumed stellar atmospheres is given in Section [6j The main focus 
of Section [7] is the comparison of fluorescence enhancement to the other major deviation to 
case B, namely collisional excitation of H°. Finally, a short summary of the paper is given 
in Section [HI 



2. The mechanism 



The abundance analysis of gaseous nebulae is based on a very simple principle. The 
intensity of a line with respect to, e.g., H/5 is given by the abundance of the emitting ion 
relative to H + times a mathematical expression describing the physical mechanisms of line 
formation involved. If such mechanisms are known and all the parameters they depend on 
(such as the electronic temperature T e ) can be determined, then the observed value of a line 
ratio can be used to measure the ionic abundance relative to hydrogen. 

The specific mechanisms involved depend on the ion and the regime considered. To 
mention but a few examples, it is known that He I lines in H II regions mainly form in the 
recombination cascade of He + , plus a density-dependent contribution by collisional excitation 
from the metastable level 2 3 S, which affects mainly the triplet lines (jOsterbrock fc Ferland 
20061 ). He II lines form in the recombination cascade of He ++ . Heavy element lines mainly 
form through collisional excitation from the ground state of abundant ionic stages. Finally, 
H I lines - a standard reference for the intensity of the other lines in the optical range - form 
mostly in the radiative cascade following recombination of H + under case B conditions. 

It is also known that case B is only a limiting case and not the exact description of 
conditions in H II regions; strictly speaking, it is not even self- consistent, as it provides no 
mechanism to depopulate the 2p level (see, e.g., iBurgessI Il958l ). However, the predictions 



[see, e.g 

made by assuming case B have generally proved to be accurate, and only in recent times 
has it become possible (and relevant) to point out exceptions to it. One such exception 
is the occurrence of collisional excitation from the ground state, which has recently been 
rediscovered and closely st udied as a consequence of the increasing need for accuracy in 



abundance determina tions (jStasiriska fc Izotovl l200ll ; iLuridiana et al.l 120031 ; iPeimbert et al. 



20071 ; iLuridianal 120081 ) . A further way in which hydrogen may deviate from case B behavior, 
which is the subject of this paper, is continuum pumping in the Lyman lines and the related 
Balmer fluorescent enhancement. 
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Continuum pumping is the excitation of a transition by means of the absorption of a 
stellar continuum photon of the appropriate wavelength. Such excitation is followed by a 
radiative cascade which, in optically thick environments such as those of H II regions, ends 
up producing a Balmer series photon plus either Lyct or 2-7 continuum. In such a way, a 
Ly/5 absorption at 1026 A can lead to the emission of an Ha photon; a Ly7 absorption at 973 
A can lead to the emission of either H/3 or Paa plus Ha; and so on. These Balmer photons 
appear independently of, and in addition to, Balmer photons formed in the recombination 
cascade of H + , so these two contributions add up in the total spectrum. As for Lya photons, 
they are scattered and gradually shifted in frequency until they escape the nebula, or are 
absorbed by dust; in dustless nebulae, a significant fraction of Lya photons is converted into 
2-7 continuum. Dust can also destroy higher Lyman photons, thus inhibiting pumping; this 
possibility will be briefly discussed in Section [7J 

It must be stressed that this Balmer enhancement mechanism is triggered by pumping 
at the Lyman-line wavelengths, i.e. it depends on the shape of the non-ionizing stellar field. 
This runs counter the common wisdom that, in a photoionized object, only the ionizing part 
of the stellar field matters; after all, the ionization structure is governed by the ionizing 
radiation, and the temperature structure is determined by the balance between heating due 
to photoionization and cooling due to collisional excitation of abundant metals and, in less 
metal-rich environments, Lya emission. As a consequence, in photoionization modeling the 
effect of non-ionizing stellar radiation on the nebula is generally neglected. Yet, the emission 
spectrum of an H II region can be altered by continuum pumping even if this mechanism 
has no other effect on the structure of the nebula. In the next section, we will justify this 
statement by d iscussing t he rol e of stellar continuum pumping in case C nebulae, referring 



in particular to iFerlandf s (119991 ) results. The following sections will be dedicated to the role 
of stellar pumping in the optically thick case, which we have labelled above "case D". In 
particular, we will estimate the contribution of stellar pumping^ on Balmer intensities and 
show that it cannot be neglected. 



2.1. Case C 

A good starting point to case D is case C, which applies to nebulae optically thin in 
the Lyman lines irradiated by a stellar continuum. In an optically thin environment, stellar 



1 A note on nomenclature: although continuum pumping and fluorescent excitation are different processes 
(the former referring to level excitation, the latter to line emission), they are linked together, so we will use 
the two terms interchangeably whenever appropriate. 
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continuum pumping results in a significant enhancement of Lya emission, sufficient to bring 
the Lyq/H/ j ratio well above the case B value, rather than below it as predicted for case 



A. iFerlandl (119991 ) explains this by showing that, in the vicinity of the star (i.e., before 
stellar Lyman-line photons have substantially been absorbed) continuum pumping provides 
a source of Balmer photons comparable to recombination. If the gas slab is optically thin, 
recombination and fluorescence act concurrently in the whole volume, so the total emissivity 
in Lya is largely enhanced by the latter mechanism. 

This argument anticipates one reason why case D has never been taken into considera- 
tion: if the gas thickness is large (i.e. the optical depth is large in the Lyman continuum), 
Lyman-line photons are used up much earlier than Lyman continuum photons, so that the 
relative effect of continuum pumping on the total emitted spectrum is diluted; the larger the 
optical depth, the more diluted the effect, until it possibly becomes negligible altogether. If 
this were always the D would not have any actual applications. There are, how- 

ever, three possible exceptions to this conclusion: a) a high accuracy in abundance analysis 
is required; b) the nebular region affected by fluorescence has an unusually large weight in 
the observation, due to aperture effects; c) the balance between the Lyman-line photons and 
the Lyman -continuum photons is more in favour of the former than in the case considered 



by iFerland . 



The first circumstance is certainly met in primordial helium abundance determinations, 
a goal that requires knowing the line formation mechanism of H and He lines with an accuracy 
better than 1%. The second circumstance is met in spatially-resolved observations sampling 
small regions in the vicinity of the ionizing cluster; observations of this kind are more and 
more common as technology improves. The third circumstan ce does n ot see m at first glance 



very probable when considering that the input spectrum of IFerlandl (119991 ) is a blackbody; 
these results, one may be tempted to conclude, overestimate rather than underestimate the 
amount of continuum pumping, since a blackbody has no absorption features. Although this 
is true, it must be con sidered that a b lackbody has no Lyman jump either, so the Lyman line- 



to-continuum flux of IFerlandl (119991 ) might actually be an underestimate. To demonstrate 
this point, in Figure [T] we compare the T e g- = 40, 000 K blackbody of that paper to the 
synthetic spectrum of a t = 3 Myr population emitting the same rate of ionizing photons 
(i.e., with the same potential to produce Ha photons by recombination; more details on this 
spectrum will be given in Section 12. 2p . The figure shows that the blackbody does not lie 
above the synthesis spectrum: rather, it crosses the near-bottom of the relevant absorption 
lines, thus underestimating the number of photons capable of continuum pumping. 



Concluding, we see that there is no a priori reason to discard the relevance of case D. 
In the following section, we will quantify in a specific case the correction represented by case 
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D with respect to case B, and discuss how this result is expected to depend on the relevant 
H II region's parameters. 



2.2. Case D 



The upper panel of Figure [2] shows the physical conditions in a model giant H II region, 
name d Ml, computed with Cloudy (version C08.00; iPorter. Ferland. van Hoof, fc Williams 

20081 ) . This versi on of the code includes a full model of the hydrogen atom with strong radia- 

tion f ield physics (IFerland fc Reeslll988l ; iFerguson fc Ferlandlll997l : IPorter. Ferland. van Hoof, fc Williams 
20081 . and Appendix [A] of this paper). For the sake of simplicity, we assumed a dustless neb- 
ula with a simple spherical geometry, with starting radius Ri n = 1 pc, hydrogen number 
density Nh = 100 cm -3 and filling factor e = 0.01. Although a blister model would be more 
appropriate to represent many, if not most, H II regions, we prefer to stick to the simplest 
possible case. The volume-integrated predictions presented in this paper do not depend on 
the assumed geometry, while the spatially-resolved predictions do; it is therefore impossible 
to give figures suitable for any geometry and any line of sight, so the reader is warned that 
our spatially-resolved predictions are given more as a rough indication of the magnitude of 
the effect than as a definitive estimate. The model's total metallicity is Z gas = 0.001, and 
the detailed elementary abundances are given in Table [21 dust is not important at such low 
metallicities. Thermal line widths are assumed; the influence of a non-thermal contribution 
is discussed in Section 15.41 The ionizing source is the spectral energy distribution (SED) of 
a t = 3 Myr stellar population, emitting ionizing photons at a rate Log Q(H°) = 52.09 s -1 . 
the SED has been computed with Sed@ (jCervino fc Luridianall2006l . and references therein), 
and does not include stellar winds. The effect of stellar winds on the predictions is discussed 
in Section O 

Assuming that the model's emission is described by case B, a fairly constant Ha emis- 
sivity profile is expected, given the shallow temperature gradient. This, however, is not the 
case: two exceptions to case B behavior can be seen in the Ha profile (Figure [2J bottom 
panel) . 

In the external part (r ^ 70 pc) there is a small bump in the emissivity profile: this 
bump is caused by collisional excitation of neutral hydrogen, a mechanism ignored in case 
B and which only becomes effective in the outer part, thanks to the drop in the hydrogen 
ionization fraction. A very low metallcity gas, which we consider here, will have a high 



kinetic temperature so these co 



this can be found in iLuridianal (120081 ) 



lisional processes are especially important; more details on 
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Fig. 1. — Comparison between the SEDs of a population synthesis model and a blackbody 
at T e g = 40 kK, normalized to the same number of ionizing photons. The position of the 
first four Lyman lines is also shown. 
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Fig. 2. — Upper panel: physical conditions in a model low-metallicity giant H II region 
ionized by a 3 Myr stellar population (see Section [2] for details). Lower panel: normalized 
Ha emissivity profile of the model, compared to case A and B predictions. 
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In the inner part there is a much larger bump, extending through a substantial fraction 
of the total radius and whose maximum value more than doubles the expected case B value, 
which is shown as a dotted line. When the upper and lower panels of Figure [2] are simulta- 
neously considered, it becomes clear that density and temperature cannot be responsible for 
the excess emissivity, as one can find pairs of points with the same T e and n e values but very 
different emissivities. In other words, the cause of this behavior cannot be sought in the gas: 
this leaves the radiation field as the only alternative. Indeed, the excess emissivity is caused 
by fluorescent excitement of Ha by continuum pumping of Lyman photons (mainly Ly/3). 
This is consistent with the position of the bump, close to the ionizing source (note that this 
bump has nothing to do with the breaks in the n e and T e g gradients seen at 6 pc, which are 
caused by the recombination of He ++ and 3+ ). 

In the following sections, we will discuss whether the behavior of this model is repre- 
sentative of real nebulae and which its consequences are for the analysis of H II regions, 
considering several specific issues. We will start by considering a simple model of a spherical 
H II region, ionized by a population synthesis model (Section [3]), and show the effect of 
continuum pumping on the Balmer intensities and ratios. In the following sections we will 
complicate this basic scenario by discussing the effect of dust, turbulence, geometry, or a 
different metallicity (Section EJ). Since the amount of continuum pumping depends on the de- 
tailed shape of the spectrum, we will also discuss the effect of assuming stellar atmospheres 
with shocked stellar winds. Since no high-resolution synthesis models are available with 
shocked winds, we will consider the case of a model H II region ionized by individual stars 
(Section EJ). Sections [7| and [S] will be devoted to the discussion and conclusions, respectively. 



3. Hydrogen fluorescent excitation predicted with population synthesis models 

The black thick curve in Figure [3] shows the near-UV portion of the SED fed into the 
model described in Section \2.2i we stress again that this SED, as all the synthesis models 
used in this work, has been computed with stellar spectra with no emission lines. In this 
wavelength range, the spectral sampling of this SED is AA = 10 A; this value is determined 
by the grid of stellar libraries used, i.e. it does not depend on the way in which population 
synthesis is carried out; any synthesis model using the same stellar l ibrary would have th e 



same resolution (this is the case, for example, of Starburst99 models: iLeitherer et al.l 119991 ). 
At such resolution, stellar absorption lines are smoothed out, so that the number of Lyman 
photons in the spectrum is overestimated. Since the mechanism studied in this work depends 
directly on the number of Lyman-line photons emitted, it is important to establish whether 
fluorescent enhancement of Ha can actually occur in H II regions or it is just an artifact 
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of low-resolution synthesis models. The only way to answer this question is by feeding a 
high-resolution SED into the photoionization code. 



3.1. Photoionization models with high resolution SEDs 



To address the issue mentioned in the previous section, we computed a new photoion- 
ization model, labelled M2, identical to Ml in all respects except for the stellar source, which 
has been replaced with a high-resolution SED in the 900-2000 A range. This SED represents 
a population of the same stellar parameters (age and metallicity) as those of the original 
SED, but has a much better spectral sampling AA = 0.3 A in the near-UV range; a portion 
of this spectrum is shown in Figure [3] (thin grey curve). The SE D is taken from a library o f 
high-resolution synthesis models computed with the code Sed@ (ICervino fc Luridianall2006l ) 
to carry out th is work. The library has been computed w i th the non-rotating trac ks of the 
Geneva group JSchaller et alill992l : ISchaerer et alill993al lbl: [ Charbonnel et al.lll993l ). includ- 
ing standard mass-loss rates computed at five metallicities: Z=0.040, 0.020, 0.008, 0.004 and 
0.001. The isochrones are calculated by means of a parabolic interpolation bet ween the tracks 



(log M , log tk) in the HR diagram on a variable mass grid as prescribed in ICervino et al. 
(120011 ) . The isochrones were generated at specific ages with a time step At = 1 Mv r for 1-10 
Myr. We have assumed a power-law IMF, with slope a = —2.35 (ISalpeterl Il955l ) and low- 
and high-mass cutoffs M\ ow = 0.1 M Q and M up = 120 M , respectively. 

To model the SED, we have combined two different sets of atmosphere models, with 
low and high resolution respectively. Low resolution atmosphere models are standard Ku- 
rucz LTE atmospheres. High resolution atmosphere models cover the wavelength range 
from 900 to 2000 A, and have been computed using the same atmospheres as the ones de- 
scribed in lMartins et al.l (120051 ); iGonzalez Delgado et al.l (120051 ) but in the UV domain. Hot 
stars (27500 < T c g < 550 00 K) are modelled as line-blanketed, non-LTE, plane-parallel, 
hydrostatic atmospheres (ILanz fc Hubeny 120031) . and s tars with 8000 < T cff < 27000 K are 
modelled with Kurucz LTE atmospheres (IKuruczl 119931 ). We have assumed no contribution 
in the UV from stars with lower tempe ratures. The high-resolut i on sy nthetic spectra are 
obtained with the program SYNTHE by lHubeny. Lanz. fc Jefferyl (119951 ). 



The portions of the SED shortward of 900 A and longward of 2000 A are equal in the 
low- and high-resolution models. The main features of these models are given in Table [3J 
together with those of the models presented in the remainder of the paper. 

A summary of the difference between Ml and M2 from the point of view of the fluorescent 
contribution to Ha is given in Table HI The three rows in the table correspond to different 
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Fig. 3. — UV spectral energy distribution of a 3 Myr, Z = 0.001 stellar population, in low 
resolution (thick black line, spectral sampling = 10 A) and high resolution (thin grey line, 
spectral sampling = 0.3 A). Also shown are the first four Lyman lines, from Lya at A = 
1215 A (which does not produce any Balmer photon) through Ly5 at 950 A. 
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ways of integrating the emissivity along the shell structure of the model: in the ID integration 
the emissivity of each shell is weighted by the thickness of the shell dr; in the 2D integration 
it is weighted by the ring surface 27rrdr; and in the 3D integration it is weighted by the 
shell volume 4nr 2 dr. Observationally, the three cases correspond to a beam-like aperture, 
a narrow slit, and an aperture including the whole object respectively. Since fluorescent 
excitation takes place in the inner part of the nebula, we expect that the effect on integrated 
luminosities of differences in the local emissivity is amplified in ID integrations with respect 
to 3D integrations. As expected, Table 0] shows that Ha is spuriously enhanced in Ml with 
respect to M2 as a consequence of the poor resolution of its input spectrum, which partially 
suppresses stellar absorption lines. The bias amounts to as much as 0.04 dex (9%) in the 
case of a beam-like simulated aperture, relative to the total Ha intensity. 
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Table 2: Elementary chemical abundances in the reference model, given as 12 + Log(X/H). 
The remaining elements have not been included in the calculation since they are quite unim- 
portant. 



Element 


Abundance 


He 


10.89 


C 


7.62 


N 


6.38 





7.87 


Ne 


7.15 


S 


6.47 


CI 


4.48 


Ar 


5.50 


Fe 


5.81 



Table 3: Main features of the photoionization models used in this work. 

Model Wavelength range (A) Nebular boundary Excitation mechanism 11 



> 2000 2000 - 900 < 900 F R CE 



Ml 


LR 


LR 


LR 


Radiation bounded + = 




M2 


LR 


HR b 


LR 


Radiation bounded = = 




M3 


LR 


suppressed c 


LR 


Radiation bounded — = 




M4 


LR 


HR fc 


LR 


Density bounded = = 




M5 


LR 


suppressed c 


LR 


Density bounded — = 





a F=fluorescence; R=recombination; CE=collisional excitation; the symbols '+', '=' and ' — ' denote that the 
mechanism is overestimated, properly included, and underestimated respectively. 
Sampled to AA/A = 0.0005 (see Section El]). 
c In the 910-1030 A range. 
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Since predicting Balmer line intensities is not a common application of photoionization 
models, which is the reason for pointing out this bias? The main reason to do it is estimating 
the amount of fluorescent excitation in real objects so that it can be taken into account 
when chemical abundances are determined. The exercise realized in this section provides 
a necessary piece of information in order to know whether such estimation can be carried 
out with low-resolution models, or high-resolution models are needed instead. The results 
above indicate that the bias in Ha introduced by low-resolution models is substantial if 
compared to our best possible guess of the total emitted Ha intensity, which is given by 
M2. Therefore, the results of Table H] indicate that low-resolution synthesis spectra are 
not suitable to study the effect of fluorescence, and that high-resolution spectra are needed 
instead; they also indicate that the predicted spatially-resolved line ratios relative to H/3 
predicted are grossly overestimated in the direction of the center; this must be taken into 
account when p hotoionization models are us ed to fit spatially-resolved spectroscopic data, 



as done, e.g., in iLuridiana fc Peimbertl (120011 ). 



Note that Cloudy automatically rebins any input spectrum, so that to take advantage 
of the larger spectral resolution it has been necessary to increase Cloudy's default continuum 
mesh, as explained in the code's documentation (Hazy, part 3, page 410). We carried out 
several different experiments in this direction exploring the resolution range 0.005 < AA/A < 
0.0004, corresponding to spectral samplings 5A > AA > 0.4A at A = 1000 A. Some of the 
spectra used internally by Cloudy after such rebinning are shown in Figure |4] in the Ly/3 
wavelength region. We found no significant differences in the emission-line spectrum below 
AA/A < 0.0005. 

As a final comment, we note that, since real massive stars may develop P-Cygni line 
profiles due to the presence of strong stellar winds, at very high spectral resolution one would 
also need to consider both the detailed line profile and the relative velocities of the star and 
the nebula, which determine which part of the spectrum fall at the relevant wavelengths. 



3.2. Net contribution of fluorescent excitation to Balmer lines 

Having established that low-resolution spectra are unsuitable to compute the fluorescent 
contribution to Balmer lines, we go on to address this problem by means of high-resolution 
spectra, considering the specific case of model M2. The fluorescent intensities in M2 can 
be estimated by computing a new model M3, identical to M2 in all respects except for 
fluorescent excitation, which is suppressed by setting to zero the relevant portion of the 
continuum (but it can also be switched off by the Cloudy command "no induced processes"). 
The fluorescent enhancement of M2 is computed as the difference between the Balmer line 
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Table 4: Bias in the total Ha intensity predicted with a low-resolution stellar continuum. 
The meaning of labels ID, 2D and 3D is explained in Section I3~T1 



Integration 


7(Ha) M i ^ 
J(Ha) M2 


ID 


0.09 


2D 


0.03 


3D 


0.02 




Fig. 4. — HR synthesis model rebinned by Cloudy to different values of the spectral resolution 
r = AA/A. 
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intensities of M2 and M3 (TableEJ)- For the reasons explained in Section lBTTl the enhancement 
is particularly large in the beam-like integration: specifically, in our models as much as 10% 
of the total intensity in Ha and 14% in Hf3 seen in the direction of the nebular center is caused 
by fluorescence rather than recombination. For completeness, we have also computed the 
fluorescent enhancement of Balmer lines of model Ml, to quantify the fraction of emissivity 
that would - erroneously - be attributed to fluorescence enhancement, had we used model Ml 
to this aim: in that case, the estimated fluorescent enhancement is almost twice as large as in 
the unbiased case (Table EJ). As anticipated in Section I3TTI this demonstrates that inputting 
low-resolution spectra to photoionization models grossly biases the estimated fluorescence 
contribution. 

Thus, the first general result is that fluorescence does indeed has a sizeable influence on 
the emitted Balmer spectrum. In the remainder of this section we will discuss the implica- 
tions of such influence. In the following sections we will briefly discuss how our results depend 
on the model's parameters, such as the geometrical assumptions on the nebular distribution 
and the parameters that define the stellar populations. 

3.2.1. Effect of fluorescent enhancement on abundance determinations 

The results listed in Table have important implications for abundance determinations. 
When a beam aperture is applied to the center of model M3, at least 10% of the total 
Balmer intensity comes from fluorescent excitation. If model M3 were a real object analyzed 
under the assumption of case B emissivity, and assuming that no fluorescent excitation is 
taking place in other ions, the inferred abundances of helium and heavy elements would be 
underestimated by up to 0.06 dex. Even in the most favourable case, the large aperture 
corresponding to a 3D integration, a 0.01 dex (2%) bias would be introduced; although this 
is negligible for most abundance measurements, it would be crucially large for high-precision 
applications such as the determination of primordial helium. Although helium lines might 
also be affected by fluorescent excitation, line ratios will in general be different from what 
predicted when neglecting fluorescent excitation. 

3.2.2. Effect of fluorescent enhancement on the Balmer decrement 

A second implication concerns extinction measurements. The upper panel of Figure [5] 
shows the local Ha/H/5 radial profile of model M2. In the innermost 20 pc and the outermost 
20 pc, this ratio differs considerably from the value expected under case B. The deviation in 
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Fig. 5. — Top panel: Local Ha/H/5 intensity ratio as a function of radius in model M2, 
compared to the case B value at the model's local temperature. Bottom panel: normalized 
Ha and B./3 emissivity profiles. 
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the external region is caused by collisional excitatio n: Ha/H/3 rises sh arply because collisional 



excitation is more effective for Ha than for H/3 (ILuridianal 120081 ) . The deviation in the 
internal region is caused by fluorescence, as can be shown by an analysis of the pattern it 
follows. Fluorescent excitation of Ha and H/3 is mainly triggered by Ly/3 and Ly7 photons, 
respectively; the optical depth in Ly/3 is larger than that in Ly7, so that Ha peaks at a 
smaller radius than H/3 (Figure [5j bottom panel). As a consequence, the Ha/H/3 ratio 
initially increases, then it decreases as Ly/3 photons become exhausted and Ly7 photons 
are still available, and eventually increases again to approach case B as both Ly/3 and Ly7 
photons are exhausted. This behavior is to be expected no matter which the initial proportion 
of Ly/3 and Ly7 photons is, as it only depends on the optical depth at the two wavelengths. 
A spectrum with a larger flux in Ly/3 might produce a Ha/H/3 ratio larger than the case B 
one, but the succession between increasing and decreasing stretches would persist. 

This interpretation of the Ha/H/3 behaviour has independently been verified by running 
models identical in all respects to model Ml except for the stellar continuum, in which small 
wavelength ranges around the Lyman lines wavelengths have been selectively suppressed. For 
example, in the model with no continuum between 1000 and 1050 A the Ha/H/3 enhancement 
disappears completely, leaving only a Ha/H/3 depression around r ~ 10 pc. For brevity, the 
detailed description of these models has been omitted. 

The difference in fluorescent enhancement of different lines is illustrated by Table [5] 
for Ha, H/3 and H7. These differences are model-specific and depend on the number of 
photons available to excite the various Lyman lines (which in turn depends on the precise 
shape of the spectrum and, in a real nebula, on line broadening and Doppler shifts), and 
on the electron distribution among different sublevels of a given level, which influences the 
downward radiative paths available and which depends on the local gas density (/-mixing 
is quite ineffective at such low densities). We cannot say whether current synthesis model 
are accurate enough for these differences to be real; but, if we take them at face value, 
the straightforward consequence of this differential enhancement is a change in the intrinsic 
Balmer decrement with respect to its case B value. 

Table M quantifies this effect for model M2, by comparing its Ha/H/3 values to those of 
model M3. The case B value for these models' temperature is 2.81; the differences between 
such value and the numbers in the table are determined by the competition among three 
factors. On one hand, fluorescent excitation enhances H/3 more than Ha (Table [5]): this 
tends to decrease the Balmer decrement. On the other hand, collisional excitation enhances 



Ha more than H/3 (ILuridianal 120081 ) . tending to increase the Balmer decrement. Finally, 
fluorescent excitation takes place in the center and is therefore favoured in ID integrations, 
whereas collisional excitation takes place in the external part of the model and is favoured 
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in 3D integrations. These dependencies explain the patterns seen in the table: in model 
M2 (with a "normal" amount of continuum pumping), fluorescence bluening dominates over 
collisional reddening in the ID aperture, while collisional reddening dominates in the 3D 
aperture; in model M3 (where continuum pumping has been suppressed), fluorescence is 
quenched, so all the ratios are larger than the case B decrement, and the difference increases 
proportionally to the weight given to external regions due to the increasing role of collisional 
excitations. The consistency of this interpretation has been verified by computing density- 
bounded versions of models M2 and M3, labelled M4 and M5. In M4 and M5, the nebula is 
truncated at a radius of 70 pc, so that the region of collisional excitation is not included in 
the computation: as a result, the Ha/H/3 ratios are either lower than (M4, density-bounded 
variant of M2) or close to (M5, density-bounded variant of M3) the case B value. 



3.2.3. Observational implications of gradients in the Balmer decrement 

We can go a bit further by recalling again that the optical depth is different for different 
Lyman-line photons, in such a way that they are selectively absorbed at different depths 
into the nebula. Note that this is true no matter which the initial proportion of Lyman- 
line photons is, so that no role is played here by the caveat of Section 13.2.21 regarding the 
accuracy of the spectrum. As by-products of the absorption of Lyman photons, different 
Balmer lines are also produced in different proportions at different depths into the nebula 
(Figure [5l bottom): as a consequence, the purely fluorescent Ha/H/3 ratio is expected to 
vary as a function of depth, being in general different from the case B value (Figure [5], top). 
Therefore, the projected Ha/H/3 ratio also depends on the particular direction considered, as 
it depends on the proportion of fluorescence-excited atoms intercepted along the particular 
line-of-sight, and therefore on the distance from the exciting stars. Of course, the ampli- 
tude and scale length of this variation depend on the scale of the fluorescence zone, which 
depends on nebular geometry and is probably very thin in most objects; specifically, it will 
be shown in Section 15.2.11 that the more plane-parallel is the gas distribution, the thinner 
is the fluorescence zone. But in regions with thicker geometries, such dependence may lie 
above the spatial resolving power and be misinterpreted as a spatial variation of interstellar 
extinction. 

There are numerous observations reporting systematic changes in the Balmer decre- 
ment j_jd_OTg_diflfen^ extinction (see, 
e.g., iMaiz-Apellaniz et al.lll998t ICannon et al.l 120021 ; iKehrig et al.l 120081 ). To investigate the 



possible role of continuum pumping in such differences, we simulate spatially resolved obser- 
vations by plotting the predicted Ha/H/3 ratio along a simulated narrow slit (Figure [6]). Note 
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that this is not the same as the 2D integration, which describes the integrated emission of 
a narrow slit. Figure [6] shows that, with fluorescent excitation properly taken into account, 
the Ha/H/3 ratio actually increases from a minimum value near the center to a maximum 
in the outskirts. With no other observations available, this pattern could be interpreted as 
a variation in true extinction and bias the dust distribution inferred from the Ha/H/3 ratio. 

We believe that, for the most part and in most cases, the observed Ha/H/3 gradient 
does indeed reflect the dust distribution; furthermore, it makes physical sense that such 
distribution is correlated to the spatial position of massive stars. However, we wish to 
call attention to the fact that a part of the gradient might also be intrinsic, arising as a 
consequence of the spatial segregation of fluorescent enhancement and collisional excitation 
of H°. Even adopting a conservative stance, the evidence presented in this section should 
be interpreted as a measure of the minimum uncertainty affecting interstellar extinction 
measurements. A rough estimate based on model M2 situates such uncertainty around 
±0.02 dex in Ha/H/3, or ± 0.05 dex in C(E(3). 



4. Influence of the spatial distribution of stars 



Assuming that the synthesis model used in this work is a good approximation to an 
actual SED, the effect of continuum pumping is inescapable. On the other hand, a natural 
objection to this conclusion is that it would cause apparent abundance gradients of apprecia- 
ble amplitude, which would h ave hardly gone unn otice d. Existing ey i dence is not conclusive: 
for example, in the works by llzotov et al.l (120061 ) and iKehrig et al.l (120081 ). radial gradients 
in the inferred abundance can neither be confirmed nor dismissed. There may be several 
causes for this. First of all, the error bar in the derived abundances is of the same order of 
the effect we are concerned with. Secondly, we are considering here a spherical geometry, 
whereas there are observational evidences of a blister geometry in some resolved H II regions; 
in such a case, the region affected by fluorescence enhancement would be compressed in a 
thin shell. Finally, the spatial distribution of stars is not always centrally concentrated as in 
our naif, highly idealized model. This comment mainly concerns giant H II regions, in which 
the scale of the system exceeds the Stromgren radius of individual stars. In such cases, we 
can expect a complex geometry, which cannot be represented by our simple model. In the 
extreme case of homogeneously distributed sources, local predictions would approach the 
volume-averaged result of a centrally concentrated source. 
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Table 5: Fraction of total intensity contributed by fluorescent excitation in models M2 and 
Ml. 
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a True fluorescent enhancement. 
b Biased fluorescent enhancement. 



Note. — Total intensities include fluorescence, recombination, and collisional excitation. 



Table 6: Dependence of predicted Ha/H/3 ratios on fluorescent and collisional excitation. 
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Fig. 6. — Solid line: Ha/H/3 projected profile along the nebular midplane. Dotted line: 
Ha/H/3 projected profile computed assuming case B emission. 
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5. Dependence of fluorescent intensity on the model's parameters 

5.1. Dependence on stellar population's parameters 

In this section we will briefly discuss how fluorescent excitation, whose effectiveness 
is a function of the Lyman-line photons available in the exciting spectrum, is affected by 
the stellar population's parameters. Only a qualitative description of the input synthesis 
models will be used to this aim. The discussion below is based on models with Z = 0.02 and 
Z = 0.001 and ages between 1 and 15 Myr. 

The outcome of a comparison between different SEDs depends on which variables are 
held fixed. In the present case the comparison must be carried out at fixed Q(H°), because 
we wish to compare the Balmer line intensity due to fluorescent excitation to the intensity 
due to recombination, which is determined by Q(B°). 

5.1.1. Age of the burst 

The fraction of Ha provided by fluorescent excitation monotonously increases as age 
increases, amounting to, e.g., 3%, 4% and 6% at t = 4, 6 and 10 Myr, respectively. Note 
that this increase is triggered by the progressive decrease in the number of ionizing photons, 
which declines with age; this is only partially offset by the increase in the absolute value of 
the equivalent width of the Lyman a absorption line taking place after 3 Myr. This can be 
seen in Figure [TJ which shows a sequence of SEDs in the age range 1-10 Myr normalized to 
the number of ionizing photons. 

5.1.2. Metallicity 

Comparison of models spanning a wide metallicity range (Z ffas ,Z* = 0.001 — 0.02) 
suggests that metallicity has no effect in very young objects, and only a slight effect in 
older objects, in the sense of larger fluorescent contributions at larger metallicities. Our 
predictions for the 3D integration of models with ages in the 1-5 Myr range are listed in 
Table El Relative contributions for other apertures might be estimated on the basis of these. 
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MA) 

Fig. 7. — SEDs around Ly/3 of a sequence of synthesis models with increasing ages, normal- 
ized to the number of ionizing photons. From the bottom up: 2, 1, 3, 4, 5, 6, 7, 8, 9, and 10 
Myr. 

Table 7: Fraction of total intensity in Ha contributed by fluorescent excitation as a function 

of metallicity and age. 

Z t (Myr) 

1 3 5 
0.001 0.02 0.02 0.03 
0.004 0.02 0.02 0.04 
0.008 0.02 0.03 0.04 
0.02 0.02 0.03 0.04 
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5.2. Dependence on gas distribution 

5.2.1. Nebular geometry 

By "nebular geometry" we exclusively refer to the distinction between thick (spherical) 
and thin (plane-parallel) models. The volume-integrated fluorescent contribution does not 
depend on geometry: roughly speaking, the relative excitation rates depend on the number of 
photons available for each process (ionizing photons in the case of recombination, Lyman-line 
photons in the case of fluorescent excitation). On the other hand, the beam- integrated flu- 
orescent contribution does depend on nebular geometry: the more plane-parallel the model, 
the closer its beam- and volume-integrated line intensities. Asymptotically, the intensities 
seen through different integrations coincide. 



5.2.2. Nebular boundedness 

Density-bounded models differ from ionization bounded models in that the outermost 
part is missing. The fraction of fluorescent excitation in these models increases proportionally 
to the extension of the missing part, until the slab becomes so thin that it falls under case 
C. 



5.3. The effect of dust 

Our simulations so far have assumed a dustless nebula. This is not a realistic assumption 
in most nebulae, and the difference may be relevant in the context of this research since 
dust grains absorb UV photons. To understand this issue, we have computed a series of 
dusty models with varying metal abundances and dust content. Our models, which span 
the metallicity range 0.004 < Z gas < 0.016, show no dependence of the relative fluorescent 
contribution on dust content (although, of course, absolute fluxes are altered). The reason 
for this behaviour is the approximate constancy of dust opacity in the relevant range, which 
roughly goes from the energy of Ly/3 to a few Rydbergs. The grains used in these models 
are Orion-type, that is relatively large; using ISM-type grains, with a distribution extending 
to smaller wavelengths, UV photons would be more effectively blocked, resulting in smaller 
and dimmer nebulae. However, the fluorescent relative contribution would stay constant. 
This is a consequence of the small relative variation of dust opacity in the relevant range, 
which goes from Ly/3 to A ~ 400 A (this r ange encloses more th an 90% of the total number 



of ionizing photons); see, e.g., Figure 19 of iBaldwin et al.l (1199ll ). which shows that the total 
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variation of the absorption coefficient in the relevant range is within a factor of two, which 
defines the maximum bias that might affect our dustless results. 



5.4. The effect of non-thermal line widths 



The line widths of actual H II regions are larger than the expected thermal width. 
The origin of this non-thermal contri bution is currently unknown, althou gh it has been 
suggested that it might be turbulence (IRozas et al.lll998l ; iRelano et al.ll2005l ). Whatever its 
cause, the effect of an increased line widt h is to make co ntinuum pumping more efficient by 
increasing the line flux available (see, e.g.. lFerlandlll992l ). Since our models so far have been 
computed under the assumption of thermal line width, we will discuss here the consequences 
of assuming increasing line widths, assuming they are caused by turbulence. 

Figure [H] illustrates the effect of turbulence on our reference model M2 at different values 
of the microturbulent velocity, which is assumed Gaussian. The fluorescent component 
increases linearly with the turbulent velocity in all the range explored; specifically, from to 
f 00 km s -1 it doubles in f D integrations and triples in 2D integrations, whereas the increase 
in 3D integrations is even larger. This gives an indication of how this paper's estimates 
should be increased to take turbulence into account. At even larger velocities, the effect 
would continue to increase. 



6. Dependence of H fluorescent excitation on the stellar atmosphere details 

In this section we will consider the case of H II regions around individual stars to 
discuss how our results depend on the details of the stellar atmosphere. In Section 12.11 
we mentioned two stellar spectral features, not present in blackbodies, that may affect the 
balance between the Lyman-line photons and the Lyman-continuum photons, and, hence, 
the relative importance of fluorescent enhancement in Balmer lines: the Lyman jump and the 
stellar absorption lines. These two features affect continuum pumping in opposite directions. 
On one hand, the ratio of non-ionizing to ionizing photons depends on the size of the Lyman 
jump: hydrogen continuum pumping has a greater relative importance in a stellar atmosphere 
with a Lyman jump than in a blackbody case. On the other, absorption lines at the Lyman 
frequencies decrease the number of continuum photons available to pump, thus quenching 
continuum pumping. Both features depend on the characteristics of the star, i.e. on its 
spectral type and luminosity class, metallicity, and strength of the stellar wind; for example, 
the size of the Lyman jump is larger in an 09 V star than in an 03 V star. 
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Fig. 8. — Relative contribution of fluorescent excitation to the intensities in Ha, H/3 and 
H7as a function of microturbulent velocity. 
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The stellar flux emitted by massive stars in the spectral range between the Lyman jump 
and 1200 A has never been studied in terms of its effect on the surrounding nebular material: 
usually, only the ionizing part of the stellar fluxes is considered to be important. However, 
the arguments of Section [2] show the importance of the non-ionizing part. Fortunately, our 
knowledge of the characteristics of the stellar flux in this spectral range has enormously 
improved in the last decades thanks to the far-UV spec troscopic observations provided by 
space telescopes (mainly Copernicus and FUSE: see, e.g.. IWalborn et al.ll2002t iPellerin et al. 
2002i ) and the development o f stellar atmosphere codes for massive stars (e.g., Hillier fc Miller 
19981 : Pauldrach et al. 2001 ). We review these characteristics here and give examples of how 
they may affect the behavior of photoionization models including continuum pumping. 

To t his aim we con s idered two sources of info r matio n: the far-UV spectral atlases com- 
piled by IWalborn et al. I d2002[) and IPellerin et al.l (120021 ). and the predictions of WM-basic 
models (IPauldrach et al.ll200ll ). Below, we describe the main stellar spectral features that 
may affect case D predictions. In particular, we are interested in the Lyman jump and those 
spectral features affecting the stellar flux emitted at the Lyman wavelengths. 

In Figure M we show the spectra of WM-basic models corresponding to four different 
spectral type-luminosity class combinations, in the wavelength range 900 - 1070 A, which 
includes the Lyman jump, the relevant H Lyman lines and other stellar features of interest. 
Only this range is represented because we are only interested in those lines that can produce 
Balmer lines through continuum pumping, i.e. Ly/3 and Lyman lines of shorter wavelengths. 
Two types of WM-basic models were computed and included in the figure (see also Table [S]): 
black solid lines represent models with a smooth wind velocity law, light dashed lines rep- 
resent models including the effect of shocks in the wind. This second type of models are 
necessary to reproduce some of the lines observed in the far-UV spectra of massive stars, 
such as O VI lines (see below). 



6.1. The Lyman jump 

This feature, located at A = 911 A, is the limit between the ionizing and the non-ionizing 
part of the stellar flux. It is produced by the presence of a sharp limit in the ionization cross 
section for the ground level of H I. Therefore, its size depends on the amount of neutral 
H I atoms in the atmosphere, which is determined by the temperature and density in the 
region of the stellar atmosphere where the stellar continuum at this wavelength is formed. 
Generally speaking, the Lyman jump will be larger at lower effective temperatures and larger 
gravities (since in these cases hydrogen is expected to be more recombined). Note, however, 
that the size of the Lyman jump may also be affected by other features of massive stars' 
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Fig. 9. — WM-basic models for stars of different spectral type and luminosity class (see [8] 
for the main models' properties). Black solid lines: models with a smooth wind velocity law; 
light dashed lines: shocked-wind models. The main stellar spectral features that may affect 
case D predictions are indicated. 
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winds (see, e.g. 



G abler et al. 


1989; 


Najarro et al. 


1996) 



Since practically no observations of the Lyman jump exist (the non-ionizing part of 
the stellar flux is generally not accesible to direct observations), we have to rely on stellar 
atmosphere models. As already mentioned, the Lyman jump is larger in stars with larger 
gravities (dwarf stars) and later spectral types (Figure [9]). 



6.2. Stellar lines 



Both observations and models show that the stellar flux at the Lyman wavelengths may 
be affected not only by the H Lyman lines, but also by resonant lines of heavy elements at 
different wavelengths that may develop P-Cygni profiles if the characteristics of the stellar 
wind are appropriate. If the P-Cygni profile is completely developed and the wind termi- 
nal velocity is large enough, the saturated blueshifted absorption may reach bluewards the 
wavelength of the closest H Lyman line (see Figure [H]) • 

As illustrative examples, we discuss here the behavior of the O VI AA 1032, 1038 doublet 
and the C ill A979 lines as a function of spectral type, luminosity class and wind character- 
istics. Figure M also shows the location of the S VI lines, which may affect other H Lyman 
lines. However, this is only given for information and its behavior is not discussed here. 

The O VI AA 1032, 1038 lines may affect the spect r al reg ion of Ly/3. The spectroscopic 
catalogues by iPellerin et al.l (120021 ) and IWalborn et al.l (120021 ) show that these lines exhibit 
a single P-Cygni wind profile in early-type O stars, both in dwarfs and supergiants. In later 
O-type spectra, two distinct blueshifted absorption features are often visible and, in the case 
of supergiants, O VI lines are still visible in Bl I stars. These lines depend only weakly 
on the temperature, luminosity, and metallicity of the star, and can still be seen even at 
metallicities as low as in the SMC. 

The presence of lines of this "superionization" species in the spectrum of OB stars 
has been commonly attributed to the effect of shock-produce d X-rays on the ionization 
distribution in hot stars winds (see, e.g.. iMacfarlane et al.lll993l for a discussion concerning 
the origin of this exotic species). In fact, only stellar atmosphere models including shocks in 
the stellar wind can produce the O VI lines. 

Figure [9] shows the synthetic UV spectra of WM-basic models with and without shocks. 
In the case of 03 stars, strong O VI P-Cygni profiles are only developed if shocks are 
included in the models. Something similar occurs in the 09.5 models where, although a 
P-Cygni profile is not developed, we can still find strong extended blueshifted absorptions. 



-33- 



In all the cases, the presence of the O VI lines totally quenches the stellar continuum flux in 
Ly/3. This only happens if the terminal velocity of the wind is above ~1700 kms"l. 



Our second example is the C ill A979 line. iPellerin et al.l (120021 ) show that, in super- 
giants, this line appears as a broad P-Cygni profile with a broad absorption component in 
spectral types from 06 to 09.7, while a broad absorpt ion feature is detectable in subtypes 
09.5 to B3 V. On the other hand, I Walborn et al.l (120021 ) showed that in SMC supergiants the 
line saturates despite the lower abundances, and does not appear weaker than in their LMC 
counterparts. Figure [9] shows the behavior of this line as predicted by WM-basic models for 
the stellar cases considered. In both 09.5 models the C ill line develops a P-cygni profile 
whose blue absorption affects the flux in Ly7. Interestingly, in the 09.5V star this line is 
inhibited by the inclusion of shocks. 



6.3. Implications of the non-ionizing stellar SED on case D predictions 

In this section we discuss how the spectral features considered above make ionized 
nebulae deviate from case B to case D. To this aim we constructed a set of spherical, dustless, 
constant density photoionization models with Cloudy. The internal radius and density of the 
models were Ri = 1 pc and % = 200 cm -3 , respectively, and no turbulent line widening was 
assumed. The SEDs predicted by WM-basic for the stellar models summarized in Table [S] 
were included as input in these Cloudy models, normalized to log<5(H°) = 52.09 s _1 . 

Figures [10] through [13] show the predicted Ha and H/3 emissivities and the Ha/H/? 
ratio for different stellar cases: similarly to Figure [HI dashed and solid lines indicate results 
obtained with or without shocks, respectively, with full inclusion of fluorescent excitation, 
while dotted lines are results obtained by switching fluorescent excitation off. The behavior 
of ratios can be easily understood by considering the spectral features of each model as a 
function of its properties. 

In 03 stars, we see that while the behavior of the Ha/H/? line ratio is very close to 
case B predictions when a smooth stellar wind velocity law is assumed, it is decreased with 
respect to case B predictions when a shocked wind stellar model is considered. In the case of 
a smooth wind law, even though the Lyman jump is small and the stellar continuum flux at 
the Ly/5 and Ly7 wavelengths is dimished, the local Ha and H/? emissivities are noticeably 
affected by stellar continuum pumping in the region close to the stars, resulting in a strong 
deviation from case B predictions. However, the Ha and Hf3 variation are similar, in such a 
way that their ratio is close to case B predictions. On the other hand, in the shocked-wind 
stellar atmosphere models, the presence of the O VI P-Cygni line (see Figure [H]) leaves no Ly/3 
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photons, and hence the fluorescent contribution to Ha is negligible. Since the fluorescent 
excitation of B./3 is similar to the previous case, the combined effect is a lower Ha/H/3 ratio 
in the inner part of the nebula as compared to case B predictions. 

A similar line of reasoning can be followed to understand the other stellar cases repre- 
sented in Figures [TU1 through [131 For example, the large Ha/H/3 line ratio predicted in the 
09.51 case when a smooth wind law is used in the stellar atmosphere model calculation is 
a consequence of the large flux in Ly/3 (which triggers a large fluorescent contribution to 
Ha) and the small flux in Ly7, due to the presence of the C III line (which quenches the 
fluorescent contribution to H/3). (Note that in this case the Lyman jump is larger.) This 
large deviation from case B vanishes when the shocked-wind stellar atmosphere model is 
considered, due to the development of a broad absorption feature in the O VI lines that 
affects Ly/3. In this case, the continuum pumping at Lj/3 and Ly7 is zero (note, however, 
that the Ha and H/? emissivities slightly deviate from case B. This is produced by pumping 
in Ly<5). 

This example illustrates how deviation from case B to case D depends on the spectral 
features appearing in the spectral region considered. Generally speaking, we can conclude 
that the characteristics of the wind has a sizeable effect on the predicted Balmer line ratios. 
Models with shocked winds have lower Ha/H/? wih respect to models with smooth wind- 
velocity laws. 



6.4. The effect of the Lyman line profile 

The discussion above has made clear that the amount of Lyman continuum pumping 
in a nebula is crucially dependent on the line profiles. These depend on both the wind 
structure and the relative motion between the stellar source and the gas, so the case should 
be considered in which the Lyman lines in the stellar continuum are in emission rather than 
absorption; in such a case, the nebula could produce more excitation than expected. We 
have taken such possibility into account by computing a sequence of models with varying 
blocking factors, whose effect is to amplify the amount of continuum pum ping by a given 
facto r, which can be either smaller or larger than one (see Appendix A in iPellegrini et al. 



2003). 



The result is illustrated by Figure [T3J which shows the ratio between the total H/3 
emission and the case B H/3 emission as a function of the blocking factor in a radiation- 
bounded spherical model ionized by a stellar continuum (solid line, left-hand scale). On the 
same figure, we have plotted the predictions for the first zone of the same models, a case 
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Fig. 10. — Local emissivities, local Balmer line ratios and projected Balmer line ratios of 
photoionization models illuminated by WM-basic models of an 03 I star, with a smooth 
wind velocity law (solid) and a shocked wind (dashed), compared to those of a model with 
no fluorescent excitation (dotted). Emissivities are in units of erg s -1 cm -2 . The distances 
to the stars are normalized to the maximum nebular radius in each case. 
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Fig. 11. — As in Figure [TUl for the case of an 03 V star. 
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Fig. 12- 



As in Figure [TUl for the case of an 09.5 I star. 




Fig. 13. 



. — As in Figure [TUl for the case of an 09.5 V star. 
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corresponding to classical case C. In either case, a blocking factor larger than 1 implies that 
pumping is amplified with respect to the predicted value. It can be seen that increasing 
values of pumping trigger H/3 intensities well above predictions, particularly under case C 
conditions. 



7. Discussion 

7.1. Comparison with collisional excitation 

In this section we highlight the differences between fluorescent and collisional excitation, 
in order to provide some useful guidelines to minimize their impact on observations aimed 
at abundance determinations. Both mechanisms provide channels to increase the intensity 
of Balmer lines, but they act with different effectiveness depending on the location within 
the nebula and the physical conditions regime. The following is a quick summary of such 
analogies and differences. 

7.1.1. Why fluorescent excitation might be more important than collisional excitation 

Along any given line of sight, fluorescent excitation takes place in the nebular region 
which lies the closest to the illuminating stars. Collisional excitation, on the contrary, takes 
place preferently in the recombination zone that marks the outer limit of the H II region. 
Such region may be lacking totally or in part if the nebula is density-bounded, whereas a 
fluorescent region always exists. A further consequence of the different position in the nebula 
is that the relative contribution of fluorescence is enhanced in slit observations. Finally, 
fluorescence appears to be happening in all kind of nebulae, with no strong preferences for 
metallicities and with a tendency to increase at increasing ages, whereas collisional excitation 
is effective only in hot, young nebulae. All in all, then, fluorescence appears to be a more 
ubiquitous process. 

7.1.2. Comparison between the uncertainties affecting either process 

The integrated value of fluorescent excitation does not depend on the geometrical dis- 
tribution of the gas; neither does the integrated value of collisional excitation, provided the 
region is radiation bounded. Therefore, the two processes are pretty much geometry-free as 
long as the whole nebula is integrated, and a detailed knowledge of geometry is not necessary. 
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Fig. 14. — H/3 to t/H/3 case B as a function of blocking factor, for a model nebula ionized by a 
stellar model with a smooth wind law. Solid line, left-hand scale: radiation-bounded model 
(case B conditions). Dotted line, right-hand scale: predictions for the inner most zone (cas e 
C conditions). Note that the ratio is, in case C, much larger than 1. See also iFerlandl (119991 ). 
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Fluorescent excitation might depend much more than collisional excitation on the dis- 
tribution of stars in the case of spatially resolved observations. This circumstance is further 
complicated by the fact that it is not necessarily dominated by the hottest stars, as are 
recombination and, probably, collisional excitation. 

Both processes also depend on the illuminating SED, but collisional excitation depends 
on it in a much weaker way than fluorescent excitation. The dependence of collisional 
excitation on the SED arises on the heating power of the ionizing photons: the harder 
the average photon, the hotter the gas, so younger clusters are more collisionally excited 
than older clusters. However, a rough knowledge of the cluster's age and its SED is usually 
sufficient to this aim. On the other hand, fluorescent excitation crucially depends on the SED, 
and specifically on the detailed profile of the relevant spectral features, which determine the 
proportion of photons available to excite the various Lyman lines. Theoretically, such shape 
depends on the stellar mix, the metallicity and deta i ls of t he atmosphere. Examples of this 



are given in Section [6] and in I Gonzalez Delgado et al.l (119971 ). where it is shown how the stellar 



Lyf3 flux may be affected by a neighbouring O VI line. This extreme dependence implies that 
moderate nebular motions might drastically change the result by Doppler-shifting the line 
photons seen by the gas. In sum, uncertainties related to the SED affect the estimation of 
fluorescent excitation, while they have a much smaller impact on the evaluation of collisional 
excitation. 

As for the Ha/H/5 ratio, it can only be increased by collisions, while it can either be 
increased or decreased by fluorescence, depending on the details of the SED and the spatial 
position considered. 



7.2. Literature on case D 

As noticed in the Introduction, not much has been written on the subject of optically 
thick nebulae illuminated by stellar radiation, which we have called "case D" . There are sev- 
eral specific reasons for this, at first surprising, fact. One of them has already been mentioned 
earlier in this section: absorption lines in stellar continua are expected to quench contin- 
uum pumping. Even if the hypothesis is made that enough photons survive, a quantitative 
analysis must by necessity rely on a detailed model spectrum; if this ingredient is lacking, 
even qualitative assessments become difficult to make (see, e.g., Section I2TT1) . Only recently 
have high-resolution stellar libraries and synthesis spectra become available, making possible 
this kind of study. From this point of view, the forthcoming publication of high-resolution 
spectral synthesis models will be a giant leap forward. 
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Fo r lines of other ions, fluorescen ce is a well-known excitation mechanism (see, e.g. 



Seaton 


1968; 


Ferland 


1992: 


Lucy 


1995) 



discussed by iBenjamin et al.l ( 120021 ). In the case of hydrogen, though, fluorescence only 
provides a small fraction of the total intensity, and until recently neither theory nor obser- 
vations were accurate enough to justify an analysis at that level of precision. As an extreme 
example, consider that at the time of publication of the Baker and collaborators' series of 
papers, Gaunt factors equal to one were still considered a reasonable approximation; and in 
1959, errors on Gaunt factors as great as 20% were still acceptable (Seaton 1959). 

A final point is brought about by improvements in technology, and in particular by the 
increasingly better spatial resolution. Since fluorescent contribution is the largest in beam- 
like observations close to the stars, the availability of spatially-resolved information might 
actually signify a worsening in accuracy. This is analogous to what happens in population 
synthesis, where the straightforward modelization of stellar populations is made difficult by 
the poor sampling associated with s mall observing; elements, in s uch a way that a probabilistic 



treatment becomes imperative (i.e.. ICervino fc Luridianal 120061 ) . 



8. Summary and conclusions 

Case D is the theoretical scenario describing optically thick nebulae irradiated by stellar 
continua. In these objects, the Balmer spectrum has a fluorescent component superposed to 
the recombination component: while the ionizing part of the stellar continuum is responsible 
for the recombination component, the non- ionizing part produces a generally small but non 
negligible fluorescent component. In the same way as collisional excitation of H°, which 
is discussed in various papers cited in this work, this mechanism represents an additive 
correction to recombination, and as such it must be taken into account in abundance analysis 
to avoid biased results. 

Case D represents a correction to the more famous case B, but, B, it is a limiting 

scenario which only asymptotically describes the behavior of real nebulae. The label "case 
D" has been given for symmetry with case C, but cases C and D are inherently different 
from cases A and B in that they are not descriptions of intrinsic gas properties, but rather 
descriptions of complex systems which depend as much on the intrinsic properties of hydrogen 
as on the "boundary conditions" represented by the stellar pumping. 

To estimate the fluorescent contribution predicted under case D it is necessary to un- 
derstand how it depends on the object's physical parameters and how it is affected by the 
observational conditions. With respect to the first point, our bottom line in this work is that 
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a proper modelization requires inputting high-resolution synthesis models into the photoion- 
ization code. Even this might be only a rough approximation as Doppler shifts and various 
line-broadening mechanisms might be acting in real nebulae, thus altering the pattern of 
photons absorption; but at least it is a physically-based approximation, and a good starting 
point for further explorations. 

Regarding the influence of observational conditions, we have shown that the relative 
fluorescent intensity observed depends on the aperture. Such dependence is straightforward 
in the case of a centrally concentrated cluster, but it becomes non trivial as more compli- 
cated spatial configurations are considered. 3D models including continuum transfer and 
photoionized by high-resolution SEDs might help in these cases, but even these could be 
insufficient to account for actual complexities, as a substantial fraction of the fluorescent 
component might come from stars of stellar type later than O. When large apertures are 
considered, geometrical features bear no influence. 

A comparison with the other major deviation from pure recombination, collisional exci- 
tation, is interesting both to clarify both mechanisms and to devise observational strategies 
to deal with them, which usually means trying to minimize them. Our study suggests that 
fluorescence might be much more important than collisional excitation in two importanty 
ways: on one hand, it may contribute a larger fraction of the intensity; on the other hand, it 
is predicted to affect virtually any H II region, whereas collisional excitation is effective only 
in hot, young, metal poor, radiation-bounded H II regions. One important difference be- 
tween the two is that while collisional excitation takes place in the outer parts of H II regions, 
fluorescent excitation takes place in the vicinity of stars: the observer willing to minimize its 
contribution on the total intensity is thus placed before a dilemma, as the former mechanism 
is minimized by spatially-resolved observations, while the latter is minimized by observations 
of the complete nebula. Finally, collisional excitations always produce a reddening in the 
spectrum, whereas fluorescence might produce both a bluening and a reddening at different 
points in the nebula. Both mechanisms must be taken into account when optical spectra are 
used to study the dust distribution of observed objects, as they can introduce small biases 
in spatially-resolved extinction measurements. 

Since this is an exploratory paper, we do not attempt to establish here an exhaustive 
calibration of fluorescent excitation of Balmer lines across the whole H II regions' taxonomy; 
rather, we aim at pointing out a problem and the tools necessary to solve it. Fluorescent 
excitation adds a small but non-negligible contribution to the observed Balmer spectra of 
H II regions and, as such, it should be studied in detail and thoroughly understood since 
neglecting it might introduce biases in abundance determinations. Such biases might be 
particularly significant when determining the abundance of primordial helium, where a high 
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accuracy is required. For such application, it is also necessary to take into account the effect 
of continuum fluorescence on the helium emission lines. 
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A. Case A, B, and C in C08 



Version C08 of Cloudy has the ability to resolve atoms and ions of the H-like iso- 
electronic sequence. For most o f its hist o ry th e code has used the com pact model of the 



Ferguson & Ferlanc 
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1997) 


neered by 


Seaton 
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1959). 



configuration are proportional to the statistical weight 2(7+1). 
probabilities, such as those given by NIST, are used. 



Term-averaged transition 



With the advent of faster processors it is now possible to resolve the nl terms within 
an n configuration. The number of states that must be considered is substantially larger. 
For instance there are 210 terms within the first 20 levels of H°, so on most processors the 
linear algebra will be roughly 2 dex slower. In the current implementation it is possible 
to adjust the number of n-configurations which are resulted into /-terms. These are then 
supplemented with several additional /-mixed n configurations to "top off" the model atom. 
Only lines produced within the nl resolved levels are reported. 

The predictions of the resolved and /-mixed models can be different at low densities 
where the real /-terms are not mixed by collisions. Continuum fluorescence, the topic of this 
paper, has different effects in these two approximations, as stressed by Ferland (1999). The 
reason is that /-averaged transition probabilities between excited n-configurations and the 
ground n=l configuration are dramatically different because only np — Is transitions occur. 
At low densities the rate of decay from an excited n to ground will be overestimated if the 
/-mixed transition probabilities are used. Another complication is that the relationships 
between the Einstein A and B coefficients are not correct if the total statistical weight of 
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the upper configuration is used. All of these problems are solved with the /-resolved model 
now used in Cl oudy. The current version of Cloudy predicts an H I spectrum in excellent 
agreement with IStorey &: Hummer! (119951 ) when the number of resolved levels is made large 
enough. 

Predictions of the fluorescent excitation spectrum can have even greater differences. 
Table [9] gives a comparison of the fluorescent-excited spectrum for an optically thin cell 
of pure hydrogen gas with unit density. It is exposed to a 50,000/T blackbody with an 
ionization parameter of 0.01. The emission would be close to Case A if there were no 
external continuum fluorescing the Lyman lines. It is in the Case C limit when continuum 
excitation is included. The first column gives the wavelength of a few L yman, Balmer, and 
Paschen lines. The second column gives the IStorey &: Hummer! (119951 ) Case A intensities 
relative to their predicted H/3. The last two columns give the Case C predictions for the 
/-resolved (current Cloudy) and /-mixed (versions before C08) cases. 
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